LOW ORDER p-MODES IN A BIPOLYTROPIC MODEL OF 

THE SUN 
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ABSTRACT 

Based on the Solar Standard Model SSM we developed a solar model in hydrostatic equilibrium using two poly- 
tropes that describes both the radiative and convective zones of the solar interior. Then we apply small periodic 
and adiabatic perturbations on this bipolytropic model in order to obtain proper frequencies and proper functions. 
The frequencies obtained are in the p-modes range of low order < I < 20 which agrees with the observational 
data, particularly with the so called five minutes solar oscillations. 

Subject headings: Solar Standard Model (SSM), Lane-Emden, Non Radial Oscillations (NRO), p-modes 



1. Introduction 

Polytropic models have largely been used in the 
study of non radial adiabatic oscillations of a 
gaseous sphere [Cowling, 1941; Kopal,1949; Scu- 
flaire,1974; Tassoul,1980]. 

We have computed the first modes of a bipoly- 
tropic model whose indices n 2 = 3.85 and ri\ = 1.5 
describes both the radiative and convective zones 
respectively of the solar interior. We used Cowl- 
ing's approximation [Cowling, 1941] which reduces 
the order of the system of differential equations 
to 2 (instead of 4). The radial part of the per- 
turbation obeys equations (1) and (2) [Ledoux- 
Walraven,1958]: 
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are proper functions, I is the degree of the spherical 
harmonic, Ti is the adiabatic exponentent equal 
to |, a is the angular frequency, N is the Brunt- 
Vaisald frequency and Li is the Lamb frequency. 
The equations (1) and (2) and the boundary condi- 
tions lead to a eigenvalue problem with eigenvalue 
er 2 and this is the problem to solve. 

2. Polytropes 

By a polytrope we understand a gas of particles 
with spherical symmetry, selfgravitating, in hydro- 
static equilibrium and with the state equation: 



P = Kp 1 = Kp 



(5) 
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K and 7 are parameters that depend only of n the 
polytropic index and the mass and radius of the 
configuration. The polytrope theory developed at 
the beginnings of the XIX century, can be used 
to know the dynamical structure of a star, that 
changes its thermodynamic state through a poly- 
tropic process in which the specific heat remains 
constant. This is a true result for several regions 
of the Sun; In Fig.l we can see that the pres- 
sure and the density are used to plot 7 = 
vs x — r/PiQ in the sophisticated SSM of Bahcall 
and Pinsonneault. Three regions clearly emerge. 
In each region the SSM output is approximated 
rather well by a straight line, indicating polytropic 



1 



behavior [Hendry, 1993]. Of these three regions the 
outermost one (71 =5/3) represents the convec- 
tive zone where heat transport is achieved by adia- 
batic convection. The inner regions constitute the 
radiative zone (72 = 1.26) where heat transport 
is achieved by electromagnetic waves. Our model 
based on this characteristic of the SSM allows to 
obtain a representation of the SSM in terms of two 
polytropes. 
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Fig. 1. — 7 at the solar interior due to SSM [Bahcall et 
al.,1995]. Regions with 7 constant are described by poly- 
tropic process. The bipolytropic model is obtained if one 
assume two values for 7; 5/3 in the convective zone, and 
1.26 in the radiative zone. 



2.1. Two polytropes within the Sun 

Following Hendry [Hendry, 1993] we use £, 9 as 
the variables in the Lane-Emden equation for 
the convective zone with index m, and 77, <j> as 
the variables for the radiative zone with index n 2 . 
The parametric polytropes for the convective and 
radiative zones are respectively 



P = Kip*, with p = Ai0(O 1-8 
P = K 2 p 126 , with p = \ 2 f(v) 3 - S5 
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The main challenge is to learn how to fit these two 
polytropes together. Since the physical quantities 
P, p and M are continuous across the interface 
(not or <fi), the variables U and V given by 



V = 
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Let us start by considering the convective zone. 
We take n\ — 1.5 thus £1 = 3.6538 and 9' = 
—0.2033. Though, this polytrope is not being used 
in the vecinity of £ = 0, it is possible to consider 
all of the solutions of the Lane-Emden equation 
for n\ = 1.5. These may be generated beginnings 
at £1 with an arbitrary starting slope and inte- 
grating inwards. Solutions with starting slopes 
less negative than 9' are of particular interest (in 
the literature, they are referred as M-solutions) 
since these are the ones which intersect the poly- 
trope that represents the radiative zone. Three 
such solutions, translated into the U, V variables, 
are shown in Fig. 2. Solutions with starting slopes 
more negative than 9' (F-solutions) do not inter- 
sect the radiative polytrope and so do not need to 
be considered here. 

Knowing 0(£) and <j>(rj) we can deduce the density 
and pressure curves p(r) and P(r). The above 
model yields a central density of p c = 1.299 x 10 5 
Kgm~ 3 and a central pressure of P c = 2.237 x 10 16 
Nm' 2 . 




Fig. 2. — Solutions of the Lane-Emden equation. The 
center of the Sun as at U = 3.0, V = and its surface at 



u = o,v 



We choose 9' n = -0.0032 as a M-Solution 



are be very useful [Chandrasekhar,1939]. 



for the convective zone. 



2.2. Characteristic Frequencies 

The characteristic Lamb frequencies L and Brunt- 
Vaisala frequency N (in notation of Mullan- 
Ulrich,1981) associated with restoring forces of 
pressure and gravity in a polytrope of order n can 
be calculated as functions of radius in terms of 
the Lane-Emden function 0. For the definition of 
9, see, for example, Chandrasekhar (1939): in a 
perfect gas, 9 equals the ratio of local temperature 
to central temperature. Thus, if we define 
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Fig. 3. — Propagation Diagram for the bp model (solid 
line) and for the polytropic index 3 model (dashed line). 



where Rq and Mq are the sun mass and the sun 
radius, and G is the gravitational constant, we can 
write a normalized Brunt- Vaisala frequency as: 



oj g (no + 1) 6 d£' 



where CC is the central condensation of the poly- 
trope, defmided as the ratio of central density to 
mean density. The parameter no is the effective 
polytropic index associated with the pulsations 
[Mullan-Ulrich,1988], which are associated with 
the adiabatic exponent T from: 



equations. Space oscillations are allowed only in 
the regions where these coefficients have opposite 
signs. The limits of the oscillations regions are 
defined by 
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In the (x,to 2 ) plane, these equations define two 
curves. In Fig. 3 we have plotted them as a solid 
line for the bipolytropic model. Have been de- 
noted p-modes and g-modes the regions of this 
plane corresponding to the conditions of position 
in the star and frequency, allowing spatial oscilla- 
tions. These regions are characterized by the pos- 
sibility of existence of progressive acoustic waves 
and progressive gravity waves respectively [Scu- 
fiairc,1974]. Thus we shall refer to these regions 
as the acoustic and the gravity regions. We have 
also plotted in the same figure the frequencies of 
the first p-modes for bipolytropic model and n = 3 
model. 

The equations (1) and (2) are very convenient for 
the analytical discussion, see for example [Unno et 
al.,1989], but for the numerical computations we 
use the more appropriate form 
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We supposed V = |. 

Similarly, the value of Lamb associated with mode 
I is given by 
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where £ is the Lane-Emden coordinate. 

3. NRO in a Bipolytropic Model 

Space oscillation properties of the solutions of 
equations (1) and (2) are related to the signs of the 
coefficients given in the second members of these 
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The regularity condition at the centre, first re- 
quires that 



ui y — Iz = 



(23) 



and second, at the surface, the cancellation of the 
lagrangian perturbation of the pressure is written 

by 
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o 



(24) 



In order to determine the solution uniquely, we 
impose the normalizing condition 
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Fig. 5. — Comparison with GOLF data for I = 2 (circle) 
and I = 4 (star). The e value is given by e- 
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at the centre. With a trial value for uj 2 we inte- 
grate equations (16) and (17), with initial condi- 
tions (23) and (25) using Runge-Kutta method, 
with a step size taken from paper of Christensen- 
Dalsgaard [Christcnscn-Dalsgaard ct al.,1994]. 
Usually this solution docs not satisfy equation 
(24) and a new integration is performed with an- 
other value of to 2 . This procedure is repeated until 
equation (24) is satisfied, using a Newton- Raphson 
method to improve the value of uj 2 . 
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Fig. 4. — Some p-modes for a polytropic models. 

Figures A and B corresponds to the n = 3 model. C, 
D, E, and F arc the same modes but in the bipolytropic 
model. The units of tj n ; are fiHz. 



The radial displacement S(r) and the pressure per- 
turbation P' are periodic space functions; the vari- 
ables v(r) and w(r) vary strongly from the cen- 
ter to the surface, then it is impossible to plot 
them directly along the axes. However, the most 

appropiatc functions £ = ±log 10 (l + |^^|) and 



C = ±log 10 (l + 



, have been plotted in Fig.. 



I GM p\ 

Their signs are chosen according to the signs of the 
variables S(r) and P'. 

4. Conclusion 

Although we do not use an atmosphere model and the 
input physics is described by polytrope structure, the 
modes obtained from the bipolytropic model are close 
to the observacional data. We show in Fig. 5 a com- 
parison with GOLF data. 
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